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Abstract 

Hydrostatically pressurized circular rings confined to two dimensions (or cylinders constrained 
to have only z-independent deformations) undergo Euler type buckling when the outside pressure 
exceeds a critical value. We perform a stability analysis of rings with arc-length dependent bend- 
ing moduli and determine how weakened bending modulus segments affect the buckling critical 
pressure. Rings with a 4-fold symmetric modulation are particularly susceptible to collapse. In 
addition we study the initial post-buckling stages of the pressurized rings to determine possible 
ring folding patterns. 
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I. INTRODUCTION 



Stiff elastic rings confined to two dimensions and subjected to in plane forces, such as 
a hydrostatic pressure, have been extensively studied in the past, mainly in the context of 
the Euler-BernouUi elastica theory and bifurcations. Related problems arise for cylinders, 
aligned with the z-axis, and constrained to have only 2;-independent deformations. The 
buckling instability of such a system under uniform pressure has been approached as an 
eigenvalue problem and issues such as existence and stability of solutions have been discussed 
in several works In addition to the problem's appeal from a purely 

academic perspective, there are physical systems where such studies find direct applicability. 
One example is carbon nanotubes, whose mechanical properties under pressure have been 
investigated with the ultimate goal of manufacturing microscopic pressure gauges Q]. 

The study of stiff tubes subject to forces finds biological applications in animal veins that 
collapse under certain loads and in the marginal band of non mammalian erythrocytes and 
thrombocytes and mammalian platelets. The marginal band is a ring-shaped microtubular 
bundle at the periphery of the erythrocyte or thrombocyte that pro vides support to the 
cell membrane and helps it withstand capillary flow stresses jo], The structural and 

mechanical properties of the marginal band have not been fully understood yet, but it has 
been claimed that non-uniformity in the spontaneous curvature might play a critical role in 
determining the flattened ellipsoidal cell shape [l^ . 

Rings under forces warrant an extensive study not only because of the above applications, 
but also because they provide some intuition about their three dimensional counterparts. 
In fact, planar pressurized thermally fluctuating rings with finite bendin g rig idity have been 
used as simplified models of fluctuating vesicles such as red blood cells 13|]. In particular. 



Leibler, Singh and Fisher (LSF) have considered thermal fluctuating closed planar chains. 



including both pressure and curvature energy terms 



14| . Different aspects of the LSF closed 



2-d random walk problem in the inflated and deflated regime have been discussed in 



the 



literature, using a number of analytical and computational techniques (see e.g. [15 



16 



18 



13, 



20|). 



In most fluctuating closed chain work, the elastic properties of the chain are assumed 
constant and uniform. Heterogeneities in the LSF model have been considered in the form 
of impurities and intercalated molecules as an additional degree of freedom in the equations 
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[see e.g. 



2l| and references therein). However, planar pressurized rings with quenched 



inhomogeneities that manifest themselves in arc-length dependent elastic properties appear 
to have received little attention in the literature. 

In this work, using a combination of analytic and computational techniques, we study 
the stability of inhomogeneous, inextensible, planar pressurized rings with finite spontaneous 
curvature at zero temperature. For large rings, the approximation of inextensibility is always 
appropriate (see below). As we shall see, spontaneous curvature becomes especially inter- 
esting when rings with regions of weakened elastic strength are considered, even in the limit 
of zero pressure. In Sec. [Ill we briefly review the stability analysis of uniform pressurized 
rings , and introduce formalism that will be used later in the paper. In Sec. Illll we con- 
sider the effects of spontaneous curvature on unpressurized rings with nonuniformities. In 
Sec. II VI we derive the critical pressure for different ring bending modulus profiles and discuss 
the aspects of the bending modulus variability that contribute most to the stability of the 
ring. The most pronounced weakening arises in rings with a 4-fold symmetric modulation 
of the bending rigidity. Last, in Sec. |Vl we explore the initial stages of the post-buckling, 
non- uniform ring behavior [26] . In the Appendices we include derivations omitted from the 
main part of the article and discuss the ring buckling problem in the context of second order 
phase transitions. 



II. BUCKLING OF UNIFORM RINGS WITH A SPONTANEOUS CURVATURE 

For completeness, in this Section we briefly review the buckling of uniform pressurized 
rings, as discussed in Refs. , and calculate the critical pressure [27|. The bending 

modulus /€ of a rod scales as the fourth power of the thickness h, k, h^; whereas the 
stretching modulus as the second power [22l]. From these facts it can be shown that, for 
rings that are sufficiently large, i.e. h <^ R, the energetic cost of bending is always smaller 
than the stretching energy [28]. Note that it is always possible to deform rods inextensibly, 
which is not the case for two-dimensional shells. In this work the rings that we consider are 
assumed inextensional. 

Upon ignoring torsion and spontaneous curvature, the bending energy of an inextensional 
circular ring of radius Vc, circular cross section and an arc-length dependent bending modulus 
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k{s) reads [22]: 
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where t 



dr 

As 



is the tangent of the ring at a point r(s) = [a;(s), parametrized by the 
arc-length s, as shown in Fig. [H The ring is assumed inextensible, so in the arc-length 
parametrization, |^| = 1. The bending energy of a ring of spontaneous curvature Cq (which 
in general is Cq 7^ ^/fc) reads: 



bend 
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ds k{s) [c{s) — cqY 



(2) 



where cis) 



m 



t(s) X 



dt{s) 
ds 



is the curvature of the ring at a point r(s) = [x{s),y{s)] 
parametrized by the arc-length s, as shown in Fig. [TJ . Here m is a unit vector perpendicular 
to the plane of the ring. The curvature c(s) is given by the derivative of the unit tangent 
vector t{s) with respect to s. The work done on the system by exerting a pressure difference 
p between the inside and the outside of the ring is: 



ri r 

x{s)dy{s) — y{s)dx{s)] — p Trr 



p 
2 



(3) 



where A is the enclosed by the ring. 

Despite the simple form assumed by the area when a Cartesian parameterization is used, 
this choice actually is not the most convenient one. In order to confine the functions de- 
scribed by r(s) = [x{s),y{s)] to a set that describes inextensional rings of perimeter 27rrc, 
we need to impose not only periodic boundary conditions r(— vrrc) = r(7rrc) but also an 
inextensibility constraint |^| = 1 for all s. This constraint is trivially satisfied if we recast 
r(s) = [x{s),y{s)] in terms of 0{s), the angle between the tangent vector and the x-axis (see 
Fig. [1]). The coordinates x{s) and y{s) then read: 



x(s) 



-nrc 
s 



COS [em 

sin m)]d^. 



(4) 



This choice simphfies Eq. (|2]), since the curvature can now be rewritten in terms of the 
tangent angle as c(s) — 



ds 



1 



-E^bend = o / k{s) 



de(s) 



(5) 



-Trrc 



The pressure now couples to an area with an interesting nonlocal dependence on the field 

o{sy. 

Am] = l J ds J d^sm[e{s)-e{0]. (6) 

— TTTc —Trrc 

Note than when k,{s) = k = const, the bending energy simplifies: 



Ehend = ^ [ ds i^^J^ ) - 2tTKCo + -KKT^cl, (7) 



2 J \ ds 

—Trrc 

SO the spontaneous curvature only contributes an additive constant in this limit. For a given 
p and k{s), the shape of the ring is determined by functional minimization of the total 
energy: 

E[^i^)] = -^bond + Epj-es (8) 

subject to the constraint 

dse*^(^) = 0, (9) 



which ensures periodicity in x{s) and y{s), and 

e{-nrc) = 9inrc) + 27f (10) 

which guarantees a complete 27r circuit of 6{s) as one moves around the ring. A last con- 
dition, 6{—7crc) = — vr/2, can be used to choose arbitrarily a solution from the group of 
degenerate solutions that are connected by rigid body rotations. 

The equation for the curvature c{s), derived from minimizing the functional of Eq. ([8]) 
plus two Lagrange multiplier terms from the real and imaginary parts of Eq. (Q, is nonlinear 
and analytically challenging (see App. For constant k, (a uniform ring) it reads: 

^(f^+"i)-P-bc = (11) 



ds^ 2^ 

where b is an integration constant. Note that the spontaneous curvature Cq does not appear 
explicitly when k is constant. The problem however can be reformulated using an ansatz 



and reduced to solving a set of transcendental and algebraic equations It is found that 
below a critical pressure equal to pc = 3k/ r^, the ring retains its circular shape. However, 
as the pressure increases above Pc, the ring suffers an elliptical deformation, and ultimately 
assumes a peanut-like shape. 

The critical pressure can also be obtained by linearizing Eq. fill I) . The quantity Ac(s) = 
c{s) — I/tc is very small for an infinitesimally deformed ring, and for constant k the zero 
pressure shape is always circular. Hence we can rewrite Eq. f|TT]) in terms of Ac(s) and and 
keep only the lowest order terms: 



A particular solution of Eq. (fT2|) is a constant proportional to the right hand side of the 
equation. However since Ac is periodic and integrates to zero, the right hand side of the 
equation must vanish. This constraint fixes the integration constant b: b = ^ — pr^. 

Eq. ( fT2l) now reduces to the equation for a classical harmonic oscillator, where the arc- 



length plays the role of time. The nontrivial solution Ac ~ sin(-v/l + + const) only 



satisfies the periodic boundary conditions provided y 1 + ^ = n, where n is an integer. 
The solution n = 1 — p = does not satisfy Eq. ([9]), and the next allowed solution (n = 2) 
gives a critical pressure Pc = Sk/t^. 

The continuous rotational symmetry of the undeformed ring is broken at and beyond the 
buckling pressure, where the ring becomes ellipsoidal. This transition can be viewed as a 
continuous, second order phase transition, with pressure as the control parameter. A short 
discussion about the transition and the order parameter is included in Appendix O 

The stability study in this section holds only for constant bending modulus and constant 
spontaneous curvature. The value of the bending modulus determines the rigidity of the 
ring and the critical buckling pressure Pc- However, provided that the bending modulus is 
constant, the spontaneous curvature does not affect the behavior of the pressurized ring, 
since the cross term cq in the quadratic bending energy density in Eq. ([5]) integrates to 
a 6'(s)-independent constant. Inextensible rings with any constant spontaneous curvature 
should have the same critical pressure and identical buckling behavior. As we shall see, the 
situation is quite different for nonuniform rings. 




(12) 
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FIG. 1: Parameterization of an inextensible ring. 0{s) is the angle between the tangent at arc-length 
s and the x-axis. 

III. SPONTANEOUS CURVATURE EFFECTS IN RINGS WITH NONUNIFOR- 
MITIES WITH ZERO PRESSURE 

As discussed at the end of Sec. [Til spontaneous curvature plays no role in the buckling 
of uniform rings. However, as is well known from everyday experience, this is no longer 
true for rings with varying bending rigidity with, for example, weakened or strengthened 
areas. Even in the absence of a pressure, if the spontaneous curvature of the ring is different 
than its natural circular curvature l/vc, then the bending stresses inherent in the ring will 
be partially released as the ring assumes a deformed non-circular shape. In this section we 
discuss the p = equilibrium shape of simple model rings with varying k, and Cq ^ 1/tc. 

Henceforth, unless otherwise stated, to simplify the formalism we will measure all lengths 
in units of the circular radius r^, which is equivalent to considering rings of total length 27r. 

A. Hinged rings 

In the absence of hydrostatic pressure, the shape of a ring with an infinitely weak link 
(i.e. a rod of spontaneous curvature cq hinged at its ends) can be determined exactly, and 
the solution can be expressed in closed form in terms of elliptic integrals. A simple desktop 
approximation is a piece of writing paper, gently folded back and attached to itself with 
flexible tape (see inset of Fig. [3](a) and Fig. [2]^a)). These solutions are closely related to the 
elastica solutions discussed in Ref. j^]- However, in our case the boundary conditions are 
fixed and the control parameter is the spontaneous curvature rather than the force applied 
on a rod. 
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FIG. 2: Parametrization of rings with non-uniformities. Inset: Hinged ring with cq ^ —1. Top: 
Intersections allowed. Bottom: No intersections allowed. (a)IIinged ring. (b)Ring with strong 
sector. 

We set the origin of the arc-length s = at the infinitely weak point and orient the ring so 
that it has bilateral symmetry across the x-axis (see Fig. [21). At the infinitely weak point the 
bending energy, and hence the torque, is zero. The ring curvature at s = 0"'' (and s = 2tt~) 
thus has to equal the spontaneous curvature Cq. For the piece of initially fiat writing paper 
mentioned above, we have cq = at this point. The continuity of the tangent angle at s = vr 
and the bilateral symmetry of the ring dictate that the ring has to cross the x axis at right 
angles at s = vr. 

These considerations translate to the following boundary conditions: 



ds 



= Co 
s=0+ 

e{TT) = 37r/2. 



(13) 



Last, the bilateral symmetry and orientation that we have chosen for the ring implies ?/(0) 
yij^), so the tangent angle 0{s) needs to satisfy: 



y{Tx) -y(0) = / sin^(s)ds = 



(14) 
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To determine the tangent angle subject to the boundary conditions of Eq. (fT3|) . we need 
to functionally minimize the ring bending energy, and implement the constraint of Eq. f|T^ 
with a Lagrange multiplier. 

For Co < 1 (preferred radius of curvature Cg^ greater than the actual one r^), after 
integrating the differential equations that are derived from the functional minimization of 
Eq. ([2]), the function 9{s) is found to be: 



Ois) 



-2 am [kX{l- s/TT),l/k] 



(15) 



where am [x, k] is the Jacobi amplitude 2J]. The constants k and A, which enter the equa- 
tions as an integrating constant and a Lagrange multiplier respectively, are determined by 
the equations: 



E 
F 



cos 



COS 



-1 



A/2 
A, 



(16) 



V2A;A/ ' 

where E [x, k] is the incomplete elliptic integral of the second kind and F [x, k] is the incom- 
plete elliptic integral of the first kind [24|- For cq ^ 1 we retrieve, as expected, a perfect 
circle, as the hinge at s = has no effect. (To recover results for Vc ^ 1, simply let cq — > cqTc 
in the above formulas.) 

The above equations also hold for cq < 0, when the spontaneous curvature is negative 
(preferred radius of curvature points outwards, which can be achieved by inverting a pos- 
itive spontaneous curvature ring). Note that for a spontaneous curvature below a certain 
threshold (approximately equal to cq = —1), the solutions of Eq. [15] will correspond to 
self-intersecting rings. 

For a ring with zero spontaneous curvature, Eqs. (|T5l) and (|T6l) simplify to: 



E(A:) 



y - 2 am [A; K(A;)(1 - s/vr), 1/A;] 
Y.{k) 



(17) 
(18) 



where Kijz) is the complete elliptic integral of the first kind. For a zero spontaneous curva- 
ture ring, the angle of the cusp at the weak point reads ■j/'q = 2 [tt — ^^(0+)] = 4sin~^ k — n, 
which equals approximately ipo ~ 0.45 237r. (See Fig. M^a) for the definition of i/jq.) 
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For Co > 1 (preferred radius of curvature Cq^ less than the actual one), Eq. (|T5ll and 
Eq.([T6]) read: 



and 



e{s) = ^ - 2 am [A;A(1 - s/n) - K{k), 1/k] ; (19) 



^-k\ = E[am{K{l/k) - k\,l/k),l/k]-E{l/k) 

TTC 

nd{Xk,l/k). (20) 



2v/F^A 



where nd(x, k) = , ^ = is the Jacobi nd function 241. 

In Fig. ©(a) we plot ring profiles, obtained by evaluating the analytic expressions Eq.( fT5|l 
and Eq. fllQI) assuming an infinitely weak point and a variety of spontaneous curvatures. As 
the curvature approaches Cq = 1 the rings tend to become perfectly circular. The curvatures 
of the black, red, green and blue ring are respectively cq = —1. (inverted ring), 0.0, 0.6 and 
1.5. The dots represent the result of a tethered chain simulation with = 601 points (along 
the lines of Sec. |V]and Ref. IJ]; for clarity only a fraction of the points is displayed). Note 
that these ring profiles were obtained assuming self-intersections are possible. However, the 
minimum energy shapes discussed here do not require self-intersections, except for cq < — 1 
or very large positive cq. In these cases the ring edges near s = will have y{s) = until an 
arc-length Sc, and the overall ring shape can then be determined by adjusting the boundary 
conditions at the origin and rescaling the ring (see inset of Fig. [2] for an illustration). 

The zero spontaneous curvature case is the simple case of a flexible rod, straight in the 
absence of torques, and hinged at its ends. For comparison with the case cq = 0, the inset 
of Fig. ([3]) (a) is a picture of ring constructed from a ribbon of paper (zero spontaneous 
curvature case) and attached at the ends with a thin segment of adhesive tape. 



B. Ring with infinitely strong sector 

As in Sec. IIII At in this section we find an exact solution to the equations derived from 
functional minimization of Eq. ([2]) for rings of zero spontaneous curvature, Cq = 0. However, 
this time we consider rings with an infinitely strengthened segment {k{s) oo) of length 
2so, over which the curvature is constrained to be exactly zero (see Fig. [2](b)). 

Upon orienting the ring so that the strengthened segment coincides with the y-axis and 
the X-axis intersects the ring in the middle (see Fig.[2](b)), we immediately see that continuity 
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FIG. 3: (a) Rings with an infinitely weak point. The spontaneous curvature co is constant along 
each ring. Curvature of the black ring: cq = —1., red: 0.0, green: 0.6 and blue: 1.5.. The dots 
represent the result of a tethered chain energy minimization simulation with = 601 points (only a 
fraction of the points is displayed). Inset: photograph of a thin ribbon of paper attached at its ends 
by flexible tape, (b) Rings of zero spontaneous curvature (cq = 0) with an infinitely strengthened 
segment. The circumference of the ring measures 27r and the length of the strengthened segment 
is 2so. Blue curve: 2so = 27r/3, green vr/6, red 7r/12 , black 0.327r. Inset: photograph of a paper 
ring with a strengthened segment equal to 16% of the total ring length. (Color available online) 
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of the tangent angle 6{s) imposes the following boundary conditions: 

eiso) = n/2 

e{TT) = 3n/2. (21) 
The continuity of y{s) and bilateral symmetry of the ring now implies: 

TT 

y{so) - y{TT) = J sin6'(s)ds = -sq. (22) 

so 

Note that y(7r) = for the configuration of Fig. [2](b). The 6{s) angle that minimizes the 
energy functional now reads: 



TT 



Cl 



^(s) = - + 2am[^(s-so),A;J, (23) 

where Ci, the curvature at sq, and the integration constant k are determined by Eq. f|2ip 
and the continuity condition of Eq. fl22|) . We can recast these two equations as: 

= (24) 

n - So 



and 

K{k) 1 



\K{k) - E{k)] . (25) 



2(l-So/vr) P 

Eq. ( |25|) can be solved numerically to determine k and subsequently ci as functions of sq. 
For So = 0, the case of a ring with infinitesimally short strong sector, the solution of Eq. fl25|) 
is /c = 0, so Cl = 1 and 9{s) =71/2 + 5, the tangent angle of an undeformed circular ring. 

In Fig. ([3D(b) we plot ring profiles, obtained by evaluating the analytic expressions Eq. fl2^ 
assuming an infinitely strong sector of varying lengths and Cq = 0. As the length of the strong 
segment approaches sq = the ring tends to become perfectly circular. The length of the 
strong segment of the red, green, black and blue ring is respectively 2so = vr/12, 7r/6, 0.32 
and 27r/3. 



C. Ring with small non-uniformities 

The boundary conditions and symmetry considerations of Sec. IIII Al and IIII Bl cannot 
be used to determine the shape of a ring at zero pressure with an arbitrary non constant 
bending modulus. We now assume that ti{s) is smoothly varying, and then consider small 
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deviations from uniformity. In this case, after functional minimization of the bending energy 
with the closure constraints (see App. [A]), we get: 

[k(s)(c(s) - Co)]"" 



CiS 



CiS 



+ Ks)(c(s)-co)]' = (26) 



where the prime stands for differentiation with respect to the arc-length. If kq is the average 
value of k(s) and we assume that both 6k, = k,{s) — kq and Ac(s) are small, we can expand 
Eq. ( l26l) to linear order in 6k, and Ac to get: 

Ko{Ac" + Ac) = (co - 1){6k" + 6K) + b (27) 

where b is an integration constant. 

From the solution of Eq. (!27|) one can derive the tangent angle that determines the 
equilibrium zero pressure shape of the ring (see App. [Bl) : 

s 

e{s) = s + ai cos s + ^2 sin s + ^^'^ ~ [ d^6K{^) + Oq. (28) 

/^o J 

— IT 

ai, a2 are constants from the homogeneous part of the solution, and 6q is an integration 
constant. The periodicity constraint Eq. fixes the values of ai, a2 to (see App. [B|) : 

TT 

«! = — / ds 6k(s) sins (29) 

VTKo J 
— n 

and 

TV 

a2 = — / ds 6k{s) coss. (30) 

— 7r 

The constant 6q can be determined if we consider that the maximum value of |c — co| 
should occur at the minimum k{s). 6k always appears multiplied by cq — 1 in this small 
deformation regime. Thus, if Cq equals 1 (the curvature of the circular undeformed ring), 
then 6{s) = s + 6q, which is the tangent angle of an undeformed circular ring. Hence, if the 
spontaneous curvature cq equals the undeformed ring curvature, then at zero pressure the 
ring remains circular, even at the presence of inhomogeneities. Inverting the sign of Cq — 1 
(by, e.g., changing the spontaneous curvature from to 2) is equivalent to inverting the sign 
of 6k, as the small curvature segments exchange place with the large curvature segments. 

In Fig. m we plot six cases of rings with cq 7^ 1 and nonuniform bending modulus. We 
focus on the particularly interesting case Cq = —1, which corresponds physically to inverting 
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a circular ring (initially with cq = 1) cut, e.g., from a yogurt cup. The ring as cut (cq = 1) 
is insensitive to variations in the dimensions along the symmetry axis {6k,{s) 7^ 0), but 
these become quite pronounced when the ring is inverted (cq = —1). In all rings of Fig. H] 
Cq = —1, and n{s) with = 1 is shown below each shape. The dotted green line is the 
result of tethered chain simulations (see Sec. |V]), the solid black line (the thickness of which 
corresponds to the magnitude of the bending modulus) is the analytical result, and the blue 
line is a circular ring of the same circumference, shown for comparison. All rings coincide 
with this circular blue ring for cq = 1. Note that for a sinusoidally modulating k,{s) with 
period 2tt, the p = configuration of the ring is circular, as seen in Fig. Hl^a). The same is 
not true for all rings with k(s) periodicity equal to 27r, e.g. the ring in Fig. Sl^d). 

As seen from the bending modulus plots, the plotted rings are far from the small 6k 
regime, when we expect the analytic result to be accurate. Nevertheless, we observe excellent 
agreement between the analytic results and the simulations, so, provided that 6k, does not 
change too abruptly, Eq. ( !28l) gives a reasonable prediction for the unpressurized ring shape. 

IV. EXPANSION IN FOURIER MODES AND CRITICAL PRESSURE OF 
NONUNIFORM RINGS 

In this section we discuss "unfrustrated" pressurized rings with spontaneous curvature 
Cq equal to the undeformed ring curvature l/r^ {cq = 1 in the units of Sec. IIIip and an 
arc-length dependent bending modulus. We determine the critical buckling pressure pc for 
various bending modulus profiles and explore the properties of k{s) that control the ring 
stability. 

A. Quadratic expansion 

As before, we describe the shape of a ring with arbitrary k{s) by the tangent angle of a 
circular reference ring, 6o{s) = const + s plus a function 66{s) which describes the deviation 
from the circular shape: 

e{s) = 9o + s + 69{s). (31) 

The constant 6q can be chosen arbitrarily and just corresponds to an overall rotation of 
the ring. All lengths are henceforth measured in units of the ring radius, so the total ring 
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FIG. 4: Zero pressure equilibrium configuration of rings with cq 7^ I/tc = 1 and non constant 
bending modulus. In all cases we consider "inverted" rings with cq = —1. The blue line represents 
a circle of the same circumference, and the distorted shapes allow optimal release of bending energy. 
All rings remain circles unaffected by the variation in k{s) when cq = 1. The dotted green line is 
the result of tethered chain simulations. To illustrate the spatial variation in k(s) the thickness of 
the rings is drawn is proportional to (0.5 + 5k(s))/15. (Color available online) 

circumference equals 27r. 

The function 66{s) can be expanded in Fourier modes: 



(32) 
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Note that since S6{s) is real, the complex Fourier coefficients satisfy 

= (33) 

The zero mode just corresponds to an overall rotation of the ring, so we set it equal to 
zero. (Upon choosing the constant in Eq. (13T|) to be Oq = n/2, we would get 0{—it) = 
— 7r/2, consistent with the parameterization in Sec. IIIII ) With this parameterization the 
inextensibility constraint is satisfied automatically. However, the modes 66n are constrained 
by the periodicity conditions of Eq. which translates to a complicated nonlinear condition 
on the Fourier modes: 

°o -AT 

E se^.sdn,...6e^j,^^.^^^^ = o. (34) 

N=l ' ni,n2,...,n]sf 

Here the summation of ni . . . , n^v extends from — oo to oo, and 5^,; is the Kroenecker delta. 

In what follows we discuss rings whose unpressurized state is circular, i.e. the spontaneous 
curvature is constant and equal to the inverse of the ring radius. As discussed in Sec. Illlt the 
ground state remains circular when cq = 1, in the absence of an external pressure. A small 
deviation from the circular, natural shape of the ring, as for example in the postbuckling 
regime for the ring when p p^, thus corresponds to small Fourier coefficients 66n- 

To second order in 69, the constraint Eq. fl34l) reads: 

oo 

^^^ + 1 E ^^nSe^^n = 0. (35) 
n=— oo 

Note that 56i is the only mode that appears at linear order, and for small deviations 5d{s) we 
expect that we can simply set 69i = 56 _i = [25]. More generally 69i is actually quadratic 
in the remaining modes, due to Eq. fl35l) . 

In the general case where the the bending modulus k(s) is an arc-length dependent 
constant and Cq = 1, the bending energy in Eq. ([2]) reads: 

^bend[^(s)] = const + dsK{s) (^^^^^ • (36) 

The spontaneous curvature and the constant part of the derivative of 6{s) cancel when 
Co = r~^, so the bending energy is quadratic in 66 with no linear terms. 

The nonuniform bending modulus k{s) can be expanded in Fourier modes as: 
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The average bending modulus of the ring will be denoted as: 1^0 = ^ fl^ K,{s)ds. 
Upon expressing the bending energy in terms of n and 66 modes we get: 



Ehend[{SOn}] 



\rn\,\l\^0,l 



(38) 



As already discussed, the mode S9o is explicitly excluded because it corresponds to an overall 
rotation of the ring. In fact, allowing /, m = would not change the value of the sum since all 
the terms containing S60 are multiplied by zero. In addition, we have excluded \m\, \l\ = 1, 
since the 66±i mode is suppressed to lowest order in the nonlinearity due to the constraint 
of Eq. ([35D. 

Upon suppressing the constant in Eq. (!36l) , Eq. (!38l) can be expressed in matrix form: 



bend 

60 



(39) 



where 66 is the infinite dimensional column vector of 66n coefficients, and the matrix ele- 
ments of -Ebend are: 

-Ebend™„ =TX mn Kn-m- (40) 

As in uniform rings, we expect that when |m| ^ 1, 59^ modes will contribute very 
little to the buckling behavior and overall shape everywhere below (and slightly above) the 
buckling transition. This observation allows us to require m < M and truncate the infinite 
dimensional 66 to a 2(M — l)-component vector, so that 

S6-M+1 



66 



602 

66 



(41) 



M 

As discussed above, the amplitudes 66o and 66±i are set to zero. We shall similarly truncate 
the bending energy expansion of Eq. ( 1371) to include only Fourier modes up to \n\ < M. This 
truncation if a good approximation for smoothly varying rings, i.e. for rings with \K,rn\ i^o 
for m > M. 
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Note that, although the 56n modes are complex, the condition of Eq. (|33|) reduces the 
complex number of degrees of freedom represented by the vector 56 by half, from 4(M — 1) 
to 2(M - 1). 

The (Hermitian) truncated matrix -Ebcnd looks like: 



M^/to, M(M .. 2Mkm-2, 

{M-lfno, .. 2{M-1)km-3, 


















9avo, 

-E-bend = 4^0, — 4^4, — 6^5, .. 

2Mkm-2 

c.c. yKn, .. 3Mkm-3 



6ki, -6k5, -9/t6, 
9ko, 



(42) 



The pressure term in the total energy {p multiplies the area, given by Eq. (^) depends on 
66{s) in a more complicated, non-local fashion. To express -Epres in terms of the 66n modes, 
we need to perform a Taylor expansion of the area in Eq. ([6]) in small S6{s), substitute with 
the Fourier expansion of 66{s) and perform the double integral in the arc-length. At first 
sight in this expansion, the 66i mode seems to appear in linear order, and therefore cannot 
be disregarded. However, the S6i and 66-i modes that appear in linear order cancel with 
sums of higher order terms as dictated by Eq. fl5Sl) . 

A more convenient approach to performing the expansion in 66n modes is to write Eq. ([6]) 
for the area A of the ring as: 



A = -Im 
2 



ds / d^e 



(43) 



If we expand the quantity e*''^*-''^ in Fourier modes as: 



(44) 



we can immediately see that the periodicity constraint exactly translates to setting fi = 0. 
Notice that, in general /„ 7^ /*„ in this representation. With the help of the /„ modes, we 
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can express the area exactly as: 



\ n>2 / 

We can then straightforwardly expand in 66^ modes and collect terms of equal order. 
We therefore see that the excess area AA = A — vr is given by: 

AA = -2n l\S9,\^ + J2\^^m\^^;^^] +0{69^). (46) 

\ m>2 / 

Because of Eq. ( I35l) . is of higher order and can be disregarded [291]. We see that the 

area, up to quadratic order, is diagonal in the 66m modes. Upon defining an area difference 

^ ^ 2 

matrix as AA such that AAmn = ~^mn'^ the total elastic energy matrix of the ring 
reads: 

^ = ^bend+pAA (47) 

At zero pressure the energy is always positive irrespective of the ring deformation, so the 
matrix -Ebend is positive definite. The circular shape of the ring which corresponds to the 
trivial configuration 59{s) = 0, is a stable energetic minimum. As the pressure increases, at 
a critical pressure Pc the lowest eigenvalue of E (or eigenvalues if the lowest energy eigenstate 
is degenerate) becomes zero, and with further increase of the pressure, negative. The circular 
shape becomes unstable at a critical pressure and the ring begins to buckle. Determining the 
critical pressure for the onset of buckling is therefore equivalent to determining the pressure 
at which the matrix E first acquires a negative eigenvalue. 

This requirement can be translated to something simpler if we consider the diagonal 
matrix whose elements are K^n = ^mn ( ) • This matrix is positive definite, and if E 
is also positive definite then the product KEK should also be positive definite. But it is 
easy to see that KEK = KE^enK — pi, where J is the identity matrix. Thus determining 
the critical pressure reduces to the easier problem of finding the smallest eigenvalue of the 
Hermitian matrix B = KE\yendK. We can now proceed to diagonalize B and determine the 
critical pressure Pc for rings with different k,{s). 



B. Bending modulus parameterization 



The full space of k{s) functions that we could explore is infinite dimensional. To simplify 
matters we identify the set of parameters that we expect to be qualitatively important in 
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explaining the trends in the value of pc for a variety of inhomogeneous rings. 

We shall compare rings with the same average bending modulus kq- Assume for sim- 
plicity that the ring is carved from a substance with uniform bulk material properties. The 
bending modulus for in-plane deformations of thin rings of a rectangular cross-section is 
then proportional to k,{s) ~ hr{sYhz{s), where hr{s) is the thickness of the ring along its 
radius and hz{s) the thickness along the perpendicular to the plane dimension 22|. The 
linear mass density is proportional to dm/ds ~ hr{s)hz{s). Assuming that the s- variation 
in the bending modulus comes through variation in the thickness hz{s), the average bending 
modulus is proportional to the total mass of the ring. Similarly, keeping hz{s) constant and 
varying hr{s), kq is proportional to ~ J^^ds(^)^. Hence, for the same total amount 
of material (average hr{s)), the maximum average bending modulus can be achieved for 
rings with constant thickness. Alternatively, more mass is required for rings with the same 
average bending modulus if we increase the amplitude of the k{s) modulations. 

In this subsection we focus primarily on rings with periodic modulations in k,{s) (rings 
with quenched random variations in the bending modulus will be discussed later). For n-fold 
periodic rings, k,{s + 2TT/n) = k,{s). As an immediate consequence of the 27r/n periodicity, 
the only modes that do not vanish in the expansion of Eq. fl37|) are integral multiples of 
n. We will consider functions k{s) with only one (global) minimum within each period 
ts = [— vr/n, Tc/n) and reflection symmetry of the ring across the line that passes through the 
minimum and the center of the ring. Without loss of generality we can choose the arc-length 
origin of the ring so that k(0) = min k,{s). It follows that K,{TT/n) = k(— vr/n) = maxK(s). 
For simplicity we define Wmin = minK(s) and Wmax = maxK(s). We now reparametrize the 
function k(s) for s G in terms of a function G{x) that parametrizes the shape of the 
defect: 

G[s/a) = —. (48) 

^max ^min 

Here, a is the width of a "defect", which we take to be a weakened ring segment similar to 
the k{s) profiles in Figs. Il](a)-|l](f). By definition G{0) = 1 and G{'k /no) = G{—7T/na) = 0. 
Upon defining G{\s\ > ir/na) = we can extend the definition of G{x) to the whole real 
axis. After averaging Eq. (l48l) over the interval tg we find: 

naq = 2n '^'"'^"J , (49) 
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TABLE I: G{x) functions describing weakened sectors and their q and values. 



G{x) 


q 




t,-[ 1,1) 


1 


3/4 


exp(— (2x)"^), m even 




1 pz-m+lN 

21/'" V m 


0(1/2 -|x|) 


1 


1 



where a normalization is provided by 

q= I G{x)dx. (50) 



Another quantity that can be used to quantify the shape of the k(s) function is its standard 
deviation relative to kq, averaged around the ring: 



cr = a/ < (k(s) - > 



Upon using Eq. fHHl) . we can recast Eq. flSTj) as: 



oo 



2> (51) 



m=l 



na 271 , , 

ivr V 

where = G'(x)^dx. 

A list of three different functions G{x) describing weakened sectors and their q and go- 
values is given in Table [H As the m exponent of the gaussian type function (second row 
of the table) increases, the function G{x) changes from a smooth, slowly varying function 
to a step like one. Increasing m in the second row thus interpolates between the smooth 
variation in the first row and the abrupt one in the third. 



C. Stability of rings with a periodic bending rigidity 

In Fig. [5] we explore the critical pressure of rings, as found by diagonalizing Eq. (j47j) for 
various periodic functions k{s) obtained by varying Wmin, a, n and G{x). In Fig. [5]^a) we use 
a sinusoidal defect profile (as shown in the first line of Table [T]). The rings displayed as insets 
depict the n = 3 case. We vary w^^ri (and through it the amplitude of the modulation) and 
calculate the critical pressure at which the ring becomes unstable. In our units, Wmin = 1 
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corresponds to a uniform ring, with pc = 3, in all cases. The case n = 1 to n = 5 are shown. 
In Fig. Mjo) we use a gaussian type defect profile, with m = 2 and a = 7r/2. In Fig. [5](c) we 
again use a gaussian type profile, but with a much narrower defect width a = 7i/8. Finally 
in (d) we plot a gaussian like function with a high enough m value (m = 12) so that the 
modulation is almost step like, and a = 7i/2. 

A striking observation, as seen in Fig. [5] (and in all other cases examined, not shown 
here), is that the n = 4 ring is particularly susceptible to collapse. As we discuss in Sec- 
tion |Vl the buckling of the ring is dominated by the 602 mode. The four equispaced defects 
couple strongly with this ellipsoidal mode of buckling, and the orientation of the elliptical 
deformation is such that the weakened points of the ring coincide with the four areas of the 
ellipse where the curvature is maximally different than the preferred spontaneous curvature 
of the ring. 

We can quantify this argument within perturbation theory. For the constant bending 
modulus, unperturbed system the matrix E is diagonal, Emq = Smp^m- Since = e_m, 
all the eigenstates are doubly degenerate. We can use the two eigenstates to construct 
two linearly independent states, one purely real and one purely imaginary, that satisfy 
69^ = 69_rn- Because of this degeneracy, we cannot straightforwardly use perturbation 
theory to predict the behavior at small 6k{s). Instead, as the unperturbed state we will 
consider the diagonal matrix plus the matrix elements at -2,2 and 2,-2 (the numbering of 
the matrix elements follows that of the 66m modes, as in Eq. ( l40l) ). That matrix can be 
straightforwardly diagonalized, and the two lowest eigenvalues now are Ai = 4(ko— p/3 — |k4|) 
and A2 = 4(^0 — p/3 + |k4|). Note that a nonzero K4 will cause Ai to go negative before A2, 
and will lower the critical buckling pressure. Upon assuming that all other k^, ttl 7^ ±4, 
are small, we use second order perturbation theory in the remaining off-diagonal matrix 
elements to obtain: 



A 



1 — 




K4I) - fc2(Ko 




+ 0(4) 



(53) 



where e "^^ = K-^/Ik^ 
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Sinusoidal modulation Gaussian modulation, wide defects 




FIG. 5: Dependence of the critical pressure on defect amplitude: we plot the dependence 

of critical pressure for various ring profiles. Here tUmin is measured in units of kq, the average 

bending rigidity of the ring, (a) G{x) = £2£(z^l±i^ (i^) G{x) = exp{—{2x)'^), a = 7r/2, (c) G{x) = 

exp(— (2x)^), a = n/S, (d) G{x) = exp(— (2x)^^), a = 7r/2. Although 3-fold symmetric rings are 

shown, the data reflects calculations for n = 1 — 5. Blue line: n = 1, red line: n = 2, green line: 

n = 3, cyan line: n = 4, black line: n = 5 (Color available online). The dashed cyan line is the 

perturbation theory prediction for n = 4. Note that the curves are not monotonic with n near 

n = 4. As n increases the curves approach the n — > oo prediction pc = 3 < 1/k dashed 

magenta line. In the insets we plot the difference dp = pc — < 1/k for n = 5 (green line) to 

n = 10 (red line). The color of the inset lines changes from green to red as n increases. The dots 

in (a) are determined from the numerical simulations discussed in Sec. |Vl 
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Eq. (153|) cannot be solved analytically, but we can iteratively approximate: 

p.^3(.„-M)-3^ i''"-:::!;:7 . (54) 

k>2 '^O -J 

This approximation is only accurate for n = 4 rings, since for n=4, the dominant k(s) 
Fourier coefficient is and the above diagonalization captures the behavior of the system 
even for relatively large modulation amplitudes. We plot Eq. ( 15^ in Fig. [5] (dashed cyan 
line), and observe excellent agreement with the numerical results for w^^^ close to 1. 

The Pc(wmin) function can be calculated analytically in the ri — oo limit of very rapid 
periodic variation in k,{s) if one considers the nonlinear differential equations for a general 
k{s) as derived by a functional minimization that leads to a generalization of Eq. (1261) (see 
App. lAl), namely 

d'[<s) (c-1))] 



ds2 



+ c/d?|.W'c(c-l) + .Wc'c] 

-p-bc = 0. (55) 



Upon expanding Eq. (!55|) to linear order in Ac = c — 1 we have a generalization of Eq. (pi 
namely, 

k{s)Ac= (b + p). (56) 



d^ [k{s)Ac] / 6 



ds^ \ /t(s) 

The integration constant b is now set equal to —p, so that the constant term that appears in 
Eq. fl56l) is zero. In the limit that 1/k,{s) varies much faster than /t(s)Ac, i.e. when ^ 1 or 
when the ring buckling is dominated by the 662 mode, we can replace 1/k in the parenthesis 
by its average: 

ir/n 



271 J n 

— 7r/n 

The dashed magenta line in Fig. E](a)-(d) is the pdn — > 00) = 3(1/k)^^. Indeed, we do 
observe that for increasing n the curves approach the limiting value pc{n oc) . In the inset 
we plot the difference 5p = Pc{n 00) — Pc{n) for = 5 — 10. The green curve denotes 
the n = 5 plot and the red the n = 10 case, and the color progresses from green to red as 
n increases. Note that the difference for each value of Wmin monotonically approaches for 
large n, as expected. 

When the functional form of G{x) allows it, 3(1/ k)^^ can be computed analytically. For 
example, for a sinusoidally modulating G{x) we have: 

3(1/k)-1 = 3v/Wmm(2Ko - tfmin)- (58) 
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FIG. 6: Left: Pc plotted versus 3(1/ k)^^ for different ring profiles and different n. In our units, 
= 1 represents a uniform ring. The discretization matrix cutoff is M = 300. The n = 3 
bending modulus profile of the ring is presented in the bottom right of the graph, the colors of each 
ring corresponding to the colors of each plot on the left. For n > 4, we observe that all profiles 
collapse on the same n-dependent curve. The dashed curves corresponding to n = 6 and n = 10 are 
derived with a matrix cutoff M = 500. For better visualization each curve is displaced by Ap = 2 
{Ap = 3 for n = 10) from the curve below. Right: Plot of versus a for n = 1 — 5. (Color online.) 
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As evident from Fig. [6l for high n the critical pressure only depends on the mean of the 
inverse bending modulus. Remarkably, plotting pc versus 3(1/k)~^ for different n in Fig. [6l 
reveals that the lines for all defect widths a and shapes G{x) to a very good approximation 
collapse on the same ra-dependent curve for n > 4. This means that the critical pressure 
is only a function of and n, regardless of the detailed shape of the modulation of 

As n increases, Pc{n oo) = 3(1/k)~-'^ as discussed earlier. The plots for narrow 
defects, large n and small 3(1/k)^^ (corresponding to deep notches in the ring) deviate from 
the expected curve. This behavior is due to the finite cutoff M that we have chosen to 
numerically compute the eigenvalues oi B = KEhendK. This cutoff limits the accuracy of 
this method for abruptly changing k{s). It implies a discretization of the arc- length, and 
when M is increased, the discretized system describes better the continuum case. For narrow 
and deep defects a finer discretization is required to represent the behavior of a continuous 
ring. Indeed, the plotted dashed lines correspond to M = 500, versus M = 300 for the 
rest of the plots in Fig. [6], and we observe that the agreement with the master 3(l//t)~^ 
dependent curve becomes better. 

Because varies too abruptly, the truncated matrix diagonalization approach does 
not allow us to examine hinged rings, like the ones of Refs. jg, Q]- However, despite the 
different regime, the qualitative result that the n = 3 periodic ring has the highest stability 
still holds for hinged rings. 

Since k{s) = kq + 6k{s), for Sk{s) <^ kq, we have 



where = is the standard derivation of bending rigidity variation, and is defined 

such that = 0. On the right of Fig. [H]we plot Pc versus a, for small cr, for n = 1 — 5. 

The curves for n = 1 and n = 2 do not agree even for cr ^ 1. However, we note that 
most of the discrepancy is explained by the difference in the defect width a. For example, 
the black and green curves {a = vr/8) collapse for small a. For small defect number n, 
when the actual defect size a/n is a considerable percentage of the ring circumference, and 
the averaged quantity is not sufficient for predicting the ring stability by itself. The 

functional form G{x) of the defect has little effect on the critical pressure, at least for small 
variances a. 




(59) 
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FIG. 7: Critical pressure of rings with quenched randomness in the bending rigidity n{s). (a) 
Critical pressure plotted versus Pc{n oo) = 3(1/ k)^^ for a large variety of random rings . Red: 
bimodal randomness, large number of weakened regions. Green: bimodal randomness , small 
number of weakened regions. Black: Smooth randomness. The five asterisks indicate the points 
which deviate the most from the prediction pc = 3{1/k)^^, denoted with the dashed blue line. 
The higher blue line is the n = 3 sinusoidal modulating k,{s) curve and the lower blue line is the 
corresponding n = 4 case. The left three columns of the plotted rings are the five realizations with 
the minimum 6p = Pc — 3(l//t)^^, plotted in ascending order, and the right three columns the five 
realizations with the maximum 6p, plotted in descending order. (Color online) 

D. Rings with a randomly modulated bending rigidity 

Finally, we discuss rings with quenched randomness in k{s). For definiteness, we consider 
first rings with bimodal stepwise randomness in the bending modulus. Thus, k,{s) can take 
only two values, Wmin and Wmax, with transitions at a number of arbitrarily chosen points 
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FIG. 8: Critical pressure of rings with quenched randomness in the bending rigidity k{s). red: 
bimodal randomness, large number of weakened regions (n^^ ~ 7 — 10). green: bimodal randomness 
, small number of weakened regions {rid — 2 — 3). The asterisks indicate the points which deviate 
the most from the prediction pc = 3(1/k)^^. (a) Scatter plot of 6p versus .^4 for rings with step 
randomness. The stars indicate the minima and maxima for the two cases. For comparison, the 
diamonds indicate the ten minima and maxima for ^5. (b) Scatter plot of Sp versus |k4| for rings 
with bimodal randomness. The stars indicate the minima and maxima for the two cases. For 
comparison, the diamonds indicate the ten minima and maxima for \k^\. (Color online) 

at positions (si, S2, S3, ... , s„^). We then pass k{s) through a low pass filter to ensure that 
Km = for m > M and our continuous elasticity formalism is valid. However, our filter 
allows for cases when k,{s) is arbitrarily small at isolated regions. In this case, the value of 
3(1/ k)^^ becomes inaccurate, skewed toward smaller values. The outlying positions of the 
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five black rings with maximum Sp in Fig. [7](a) is due to such effects. 

We also generate rings with more continuous quenched randomness by discretizing the 
arc-length to ~ 600 points, and assigning a random (chosen from a box distribution) 

to each of them (we control the average spread of 6k{s) in the box distribution so that we 
can vary a). We again pass k,{s) through a low pass filter to generate smoother versions of 
these more continuously random rings. We then determine the critical pressure for particular 
realizations by matrix diagonalization, as we did for periodic rings. In Fig. El^a), the red 
data points correspond to the bimodal randomness with a large number of weakened regions 
[ud — 7 — 10), the green data to step randomness with small na — 2 — 3, and the black to local 
randomness. The critical pressure for various realizations and values of 3(l//t)^^ is plotted 
in Fig. [7](a). Note that the majority of the results lie below the Pc{n — > oo) curve, denoted 
with the dashed blue line. For fixed average inverse bending modulus periodicity 
and large n in the bending modulus increases the ring stability. The lower blue curve is 
the n = 4 periodic case, the upper is the n = 3 case. We observe that the n = 4 periodic 
case provides a lower limit for the critical pressure of rings with quenched randomness. To 
visualize the k(s) profile of rings with quenched randomness in Fig. El^b), we plot the five 
cases with the maximum Sp = Pc — 3(1/ k)^^ in descending order (right three columns), 
and the five cases with the minimum 6p in ascending order (three left columns). The rings 
with continuous randomness with the maximum 6p (red color, right column) include regions 
where the bending modulus nearly vanishes and cannot adequately represent the continuous 
behavior. 

At first inspection, the ring bending modulus profile does not seem to explain the low 
stability of the left column for bimodal randomness (red and green columns). Since the 
limiting stability values are determined by the periodic cases, we need to describe how 
close is a ring to an n-fold periodic profile. For that purpose, we define the periodicity self 
correlation 

e„ = --— , J2i^Kis)SK [s + ^p (60) 

^ ' m=l ^ ^ 

This quantity is a measure of how close is the shape of a ring to being n-periodic. For rings 
with 6k,{s) = 6k,{s + 27r/n), it reaches its maximum value ^„ = 1. As we can see in the 
scatter plot Fig. [8](a), most of the minimum Sp values lie on the bottom right part of the 
plot. Scatter plots for any other n do not exhibit such a trend. For comparison, we plot the 
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minimum and maximum 6p values plotted against ^5 as diamonds. We should note here, 
however, that the measure ^4 will be very high for an n = 8 periodic-like ring. To distinguish 
between the two cases we can look at the K4 Fourier mode as a measure of 4-periodicity. In 
Fig. [H](b) we present the scatter plot of 5p versus and we observe a similar pattern as in 
Fig.EKa). 



V. RING DEFORMATIONS BEYOND THE BUCKLING TRANSITION 

The quadratic approximation used in Sec. [IV] is sufficient for determining the stability 
of the rings, but cannot provide information about the post-buckling behavior of the rings, 
as the higher order 59 modes are necessary to describe the behavior of the system. In this 
section we briefly discuss the initial stages of the post-buckling behavior of rings with a 
periodically modulating bending modulus, studied using computer simulations. 

For the simulations we used a tethered chain model very similar to the LSF model of 



Ref. 



14| . In our case, we considered rings with spontaneous curvature cq = 1/rc so that the 



p = shape is stress free and circular. The discretized bending energy we simulate reads: 

1 ^ 

^bend = ^ 5Z S(l - COS(A0,- - A^o)) (61) 

where A6'j is the angle between the links j and j — I- K,j is the arc-length dependent 
bending rigidity characterizing the region between Sj and Sj_i. The spontaneous curvature 
enters Eq. f l^T]) through A^q which is the zero bending energy angle between two successive 
links. For cq = l/r^, A^q = ^n/N. Here, As is the tether length, and when N is large. 
As ~ 2TTrc/N. The pressure term is as in Eq. ([3]), only now the area is calculated as the 
area of an A^-sided polygon. For computational convenience, we did not strictly impose 
inextensibility to the tethers, but instead allowed them to stretch while imposing a very 
high energetic cost to deviations of the link length from the desired length As: 

N 

E,tr = Y,tr $^(|r,+i - r,| - As)^ (62) 
i=i 

with the stretching modulus chosen such that K,tr ^ ^^^^r- 

Using a simple adaptive step gradient descent method we determined the minimum energy 
ring configuration for different sinusoidally modulated bending moduli choices {kj} and 
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pressures p. At each pressure p, we decomposed the minimum energy ring profile into S6 
modes. We plot some representative cases in Fig. [91 At the left column, we plot the pre 
and post buckling shape for n = 1 — 4 and n = 10 rings, accompanied by the amplitude 
of the relevant 66^ modes as pressure increases, on the right. The post buckling shape is 
shown for the maximum pressure displayed on the right. In all cases plotted Aw = 0.8. The 
66m modes are zero in each case below a certain critical pressure, and acquire a non zero 
value above that. Due to finite size effects, the second order phase transition is rounded out 
slightly, and determination of the critical pressure from the graphs becomes harder. Critical 
slowing down makes the simulations converge slowly near the critical pressure, thus making 
the data in this region less reliable. Hence, to estimate the value of Pc from the data, we 
fitted a 2nd order polynomial to the post-buckling data 1 1 and extrapolated the fit to 
\662\ = 0. The critical values extracted this way are shown, for several Aw, in Fig. [5](a), 
and were discussed in Section IIVCI We observe excellent agreement between the numerical 
data and the eigenvalue method. 

The approximation that the 66i ring amplitude can be neglected is still valid after the 
buckling transition. 66i is always much smaller than the other modes that dominate the 
buckling, especially close to the critical point. In particular, regardless of the bending 
modulus profile, the S62 mode always dominates, driving the transition. The rings we 
have studied always buckle into an ellipsoidal shape. As an example of non-sinusoidally 
modulating ring, in Fig. [10] we plot the initial post-buckling shapes of rings with a gaussian 
type, {m = 8, a = 7r/3) bending modulus, and observe that again the 662 mode dominates. 
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FIG. 9: Left: Rings with periodic sinusoidal bending rigidity profiles above and below the critical 
buckling pressure. On the right of each ring pair, we plot the lowest nonzero 69m amplitude versus 
the applied hydrostatic pressure p. Blue line: \69i\, green line: \S92\, red line: \693\, cyan line: 
\594^\, magenta line: \S9^\. Note the lower buckling threshold for m = 4. (Color online) 

APPENDIX A: DERIVATION OF THE CURVATURE EQUATION FOR PRES- 
SURIZED RINGS 



For completeness, in this Appendix we generalize the derivation found in Ref. [3| for 
T^const and derive the non-linear differential equation for the curvature of the hydro- 
statically pressurized ring with an arc-length dependent bending modulus. 
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FIG. 10: Post buckling states of rings with abrupt step-like modulation and a = vr/S. The 
ellipsoidal mode of deformation dominates. 

The tangent angle 6{s) of the ring that minimizes the sum of the bending energy and the 
pressure term subject to the ring closure constraint is found by functional minimization of: 



Etm] 



ds k(s) 



d9{s) 



1 2 



Co 



TV S 



+ Ai / ds cos6'(s) + A2 / dssin6'(s) 



(Al) 



We have set the radius of the ring Tc = 1. The first term is the bending energy -Ebend, the 
second term is the work done by the pressure -Epres and the last two terms are the Lagrange 
multiplier constraints that come from the real and imaginary parts of Eq. iQ. The functional 



derivative 



SEt[e{s)] 
s[e{3)] 



of Eq. (lAip if zero for the energy minimizing 6{s): 



6[e{s)] 



[k{s) (e'is)~co)]' + p f d^cosieis)-eiO) 

J —n 

cos{e{s) - ^(0) - Al sin(^(s)) + A2 cos(^(s)) = (A2) 



where the prime stands for differentiation with respect to the arc-length. To eliminate the 



constants Ai and A2 we apply the operator C = qtj^-^ 
After some algebra, it is easy to see that the sum C 



9'{s) ds 



to both sides of Eq. (jM 



5Et[e{s)] 

5[e{s)] 



+ 



5Et[e(s)] 

5[e{s)] 



is independent of Ai 
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and A2 and equal to: 



5Et[e{s)] 
1 



c{s) 



ms)] 

'k{s){c{s) - Co)]' 



c{s) 



- ['^(s)(c(s) - Co)]' 



p 



C[S 



C[S 



(A3) 



where cis] 



de 

ds ' 



Eq. (|TT1) can be easily derived from the general Eq. ( ]A3[) . by setting k{s) = k =const, 
multiplying the equation by c(s) and integrating once with respect to s. The constant b in 
Eq. f fTTj) comes from the integration. Note that since in Eq. flA3P the spontaneous curvature 
Co appears only coupled to derivatives of k{s), it drops out of the equation when the bending 
modulus is a constant. 

Similarly, Eq. (1551) comes from Eq. (1A3P by setting Co = 1, multiplying by c(s) and 
integrating with respect to s. Eq. ( I26l) follows immediately from Eq. ( ]A3[) if we set p = 0. 



APPENDIX B: TANGENT ANGLE OF THE EQUILIBRIUM ZERO PRESSURE 
SHAPE 

A particular solution of the non-homogeneous Eq. (I27|) can be easily determined by 
inspection to be Aci = (co — 1)6k/k,o + B/kq. Since Aci = — 1 needs to integrate to 
zero to be a vahd solution if 9{s) is periodic, we find that 6 = 0, so Aci(s) oc 6k{s). Adding 
the solution to the homogeneous equation, the curvature of the ring finally reads: 

1 



cis) = —a\ sm s + 0^2 cos s + 



Co 



Ko 



-5k + 1, 



(Bl) 



where a\ and 02 are s-independent constants. The tangent angle d{s) of Eq. ( l28l) can be 
easily evaluated by integrating Eq. fIBll) with respect to s. 6{s) will not satisfy the closure 
constraint Eq. (Q for arbitrary constants ai and 02 • To determine these constants we 
expand Eq. iQ to lowest order in 6{s) — (s + ^0) 



«! cos s + ^2 sin s + 



(co 



e'Ms ~ 0, 



(B2) 



and solve the real and imaginary part of the equation for ai and ^2 respectively. 
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APPENDIX C: BUCKLING OF RINGS AND SECOND ORDER PHASE TRAN- 
SITIONS 



In this Appendix we discuss how the buckhng of pressurized rings can be viewed in the 
context of the Landau phenomenological theory of phase transitions. 

It is known that the buckling of a uniform, hydrostatically pressurized, inextensible ring is 
a subcritical pitchfork bifurcation 4, l5|. The stable, circular state becomes unstable above 
the critical pressure Pc, where the ellipsoidal solution branch emerges in the bifurcation 
diagram. The ellipsoidal solution continuously becomes the circular solution as p —>■ p^, 
a behavior reminiscent of a second order phase transition, with pressure being the control 
parameter. 

The energy of the ring, when k(s) = const. = kq and 9{s) = s + S9{s), reads 



2 



— TT 
TT S 

+lJ^^J d^sin [s - e + 59(3) - 59(0]- (CI) 

— TV —IT 

This quantity is similar to the free energy of a one dimensional spin model m{s) = 
[cos 6 9 (s), sin 6 9 (s)], where the spins are arranged on a line with periodic boundary con- 
ditions. However, we should proceed with this analogy with caution. The interaction of the 
spins m{s) in the second term is non-local, with the spins at positions s and s + tt being 
coupled. Another comphcation is that the ring closure constraint Eq. should be included 
in the free energy Eq. ( ICll) with two Lagrange multipliers, one for the real and one for the 
imaginary part. The modified free energy is Eq. (jAl]) with 59{s) = 9{s) -s- const. To de- 
termine the value of the Lagrange multipliers Ai and A2 and evaluate the free energy density 
functional, one could solve Eq. (1A3P and substitute the solution into Eq. (lA2p . but this is 
not the approach we will follow here. Instead, we will use the expansion of S9{s) in Fourier 
modes, discussed in Sec. HVl and implement the closure constraint through discarding the 
69i Fourier mode of the expansion. 

We first discuss the important lengthscales and parameters that control the physics of 
the problem. In this analysis, the temperature T of the system is zero. There are no 
thermal fluctuations or entropy terms and the configuration of the ring always is at a local 
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A6 modes of uniform pressurized ring 




3.15 



FIG. 11: Solid blue lines: Expansion of the analytical solution S9(s) of a pressurized uniform ring 
in A6n = 2\69n\ modes. Dashed red line: Prediction of a model where only the soft mode 502 has 
been taken into acount. Note the excellent agreement even far from the critical pressure p = 3. 



energetic minimum. The thermal persistence length i ~ KQ/ksT is infinite, the definition of 
a correlation function in the absence of temperature becomes problematic, and there is no 
obvious diverging lengthscale. However, from the bending modulus kq and the pressure p 
we can construct a "bending" lengthscale ib{p) = {f^o/pY^^ that is of relevance to this zero 
temperature system. This "bending length" , divided by the ring radius can be thought of 
as a measure of the importance of the bending term of the energy compared to the pressure 
term. For p close to pc, h{p) ~ rc, but for high pressures h{p) ^ ^c, indicating that the 
ring can have equilibrium states that are sharply bent jscj]. Note that for high t emp eratures, 
where i{T) <^ ih, the (self-avoiding) ring will behave like a branched polymer 13|, Il4|. 

Only the longest wavelength (lowest bending energy) Fourier modes of 66 (s) take part 
near the transition. This can be seen by linearizing Eq. ffTTl) close to the transition point, 
where Ac is small, or equivalently expanding the energy Eq. ICll to quadratic order. At 
p = Pc, as discussed in Sec. HV] the 602 mode becomes a soft mode, and is the one driving 
the buckling transition. 

Having established that the soft mode S62 is the one driving the transition, we can expand 
E[6d{s)] in Fourrier modes and discard all modes except 5^2- 



^ ^ 8 



TT 



^o-l) m' + ^16021' + ... ]+pn 



(C2) 



where we have set Tc = 1 for simplicity. We see that there are no third order terms coming 
from the expansion of the area. This is something that we expect, since the energy should 
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be invariant under rotations of the ring and it cannot depend on the phase of the mode 5^2- 
Eq. flC2|) for p < Skq has a minimum at \662\ = 0, whereas for p > Pc the 662 = solution 
becomes unstable. For p — Pc Pc the position of the minimum moves to 

ise,\^,ll('^)"\ (C3) 



Pc 

The exponent of Eq. (1C3P could have been easily obtained if 662 was recognized as the 
order parameter of the system, and then we constructed the "mean field" energy Eq. (lC2p 
by expanding in 662 and recognizing the symmetries of the system. In this sense, we can 
describe the buckling transition of the ring as a second order phase transition, with the 
pressure as the control parameter. In Fig. [TT]we show the expansion of the exact solution in 
Fourier modes (solid blue lines) and the prediction of Eq. flC2l) (dashed red lines). Note that 
the neglected 66^ and 66q modes are very small and we observe excellent agreement 31 1. 
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